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ABSTRACT: 

Cigarette  smoking,  the  major  cause  of  lung  cancer,  creates  a  “field  of  injury”  throughout  the  respiratory  tract.  We  have 
previously  shown  that  gene  expression  from  bronchial  epithelial  cells  reflects  the  physiologic  response  to  cigarette  smoke 
exposure  and  can  serve  as  a  diagnostic  biomarker  for  lung  cancer.  The  purpose  of  this  Idea  Development  Award  is  to 
conduct  single  cell  RNA  sequencing  on  airway  epithelial  cells  obtained  from  smokers  with  and  without  lung  cancer  to 
identify  cell-type  dependent  gene  expression  alterations  in  the  lung  cancer  field  of  injury. 

Cells  are  being  collected  by  brushing  the  right  mainstem  bronchus  of  smokers  undergoing  bronchoscopy  for  the  suspicion 
of  lung  cancer.  A  protocol  has  been  developed  to  isolate  single  cells  from  these  bronchial  brushings  using  fluorescence- 
activated  cell  sorting  (FACS).  Next,  we  have  implemented  an  adapted  version  of  the  CEL-Seq  RNA  library  preparation 
protocol  that  includes  plate-,  well-,  and  transcript-specific  barcodes  allowing  hundreds  of  cells  to  be  pooled  together  and 
sequenced.  We  have  also  developed  a  computational  pipeline  to  process  the  sequencing  data  into  gene  level  counts  for 
each  cell.  In  order  to  demonstrate  that  the  methodologies  described  are  working,  a  pilot  experiment  was  conducted  in 
which  bronchial  epithelial  cells  from  1  former  smoker  and  1  current  smoker  were  profiled  (n=24  cells  per  donor).  Data 
generated  from  this  experiment  illustrated  that  there  is  low  cell-to-cell  technical  variation  and  known  smoking-associated 
gene  expression  alterations  can  be  detected.  The  success  of  this  initial  experiment  is  significant  because  it  demonstrates 
that  all  the  protocols  are  currently  in  place  to  begin  processing  samples  collected  from  smokers  with  and  without  lung 
cancer. 


Over  the  next  year  we  plan  to  sequence  hundreds  of  cells  per  donor  from  about  30  smokers  with  and  without  lung  cancer 
to  discover  known  and  novel  cell  types  with  gene  expression  changes  associated  with  the  presence  of  lung  cancer.  These 
discoveries  may  enhance  current  lung  cancer  diagnostics  as  well  as  suggest  potential  new  therapeutics  for  lung  cancer. 


INTRODUCTION: 

Cigarette  smoking,  the  major  cause  of  lung  cancer,  creates  a  “field  of  injury”  throughout  the  respiratory  tract  by  inducing 
molecular  alterations  such  as  allelic  loss,  p53  mutations,  changes  in  promoter  methylation  and  telomerase  activity^’^.  We 
have  previously  shown  that  gene  expression  from  bronchial  epithelial  cells  reflects  the  physiologic  response  to  cigarette 
smoke  exposure®’^.  Importantly,  we  have  extended  this  airway  field  of  injury  to  the  study  of  lung  cancer,  and  have 
identified  a  bronchial  airway  gene  expression  signature  that  can  serve  as  a  diagnostic  biomarker  for  lung  cancer®  that 
performs  independently  of  clinical  risk  factors  for  disease®.  The  lung  cancer  diagnostic  biomarker  has  been  subsequently 
validated  in  in  a  large  clinical  trial^®  and  has  been  commercialized  by  Veracyte,  Inc.  and  is  known  as  PERCEPTA™. 
Advances  in  technology  for  amplification  of  low  amounts  of  RNA  combined  with  next-generation  sequencing  have 
produced  the  ability  to  characterize  the  transcriptome  of  individual  cells.  While  the  bronchial  brushings  examined  in  our 
previous  studies  have  captured  a  relatively  pure  population  of  bronchial  epithelial  cells,  we  are  unable  to  discern  which 
airway  cell  type  or  types  are  responsible  for  the  gene  expression  changes  observed  nor  characterize  gene  expression 
variation  between  cells.  Variation  in  gene  expression  across  single  cells  can  be  used  to  define  unique  subpopulations  of 
cells  that  may  be  independent  of  known  markers  or  cell  morphology  that  may  associate  with  lung  cancer.  We  hypothesize 
that  the  lung  cancer-specific  gene  expression  in  the  bronchial  epithelium  might  be  restricted  to  specific  known  cell  types 
(e.g.  basal  cells)  or  molecularly  defined  subpopulations  of  cells.  The  goal  of  this  study  will  be  to  use  single-cell  RNA 
sequencing  to  identify  cell-type  dependent  gene  expression  alterations  in  the  lung  cancer  field  of  injury  and  to  molecularly 
identify  novel  subpopulations  of  cells  that  are  associated  with  lung  cancer.  These  novel  molecular  insights  hold  the 
potential  to  improve  the  diagnostic  utility  of  the  airway  epithelium  for  lung  cancer  and  to  guide  new  therapeutic  strategies 
for  lung  cancer  prevention. 


KEYWORDS:  single  cell,  bronchial  epithelium,  mRNA  sequencing,  and  lung  cancer 


ACCOMPLISHMENTS: 

What  were  the  major  goals  of  the  project? 

•  Specific  Aim  1 :  Identify  which  cell  types  in  the  airway  epithelium  harbor  the  lung  cancer-specific  alterations  in 
biomarker  genes  by  single  cell  RNA  sequencing 

0  Major  Task  1 :  Isolate  and  sequence  the  RNA  of  single  epithelial  cells  from  the  bronchus  of  smokers  with 
and  without  lung  cancer  (n=15  subjects/group,  n=960  single  cells/subject). 

■  Subtask  1 :  Approval  of  IRB  and  HRPO  (1  -2) 

The  IRB  at  Boston  University  School  of  Medicine  approved  the  study  protocol  on  September  1 0, 
2014  but  notification  of  the  outcome  was  received  on  December  4,  2014.  The  HRPO  approval 
was  obtained  on  December  19,  2014.  This  process  took  approximately  6  months  to  complete 
and  delayed  sample  collection  by  about  4  months.  Completion  percent:  1 00% 
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■  Subtask  2:  Collection  of  airway  brushings  from  30  subjects  at  BUMC  (2-12) 

We  have  collected  airway  brushings  from  9  subjects  at  BUMC  In  the  past  6  months  after  full  study 
approval.  We  have  collected  9  brushings  from  current  and  former  smokers  under  clinical 
suspicion  of  lung  cancer.  Of  those  subjects,  5  have  been  confirmed  to  have  lung  cancer,  3  do  not 
have  lung  cancer,  and  1  has  yet  to  be  determined  (Table  1).  Additionally,  from  an  external  cohort 
(see  “Sample  Collection’’  section  under  “What  was  accomplished  under  these  goals?’’  we  have 
collected  1 1  brushings  from  healthy  current  smokers  and  1 0  brushings  from  healthy  former 
smokers. 

As  collection  has  been  slower  than  expected  we  are  In  the  process  of  collaborating  with  Dr. 

Robert  Browning,  an  Interventional  pulmonologist  at  Walter  Reed  Medical  Center,  to  collect 
additional  brushings.  We  are  currently  adding  his  center  to  our  IRB  protocol.  Completion 
Percent:  33% 


■  Subtask  3:  Sorting  of  cells  from  brushings  using  FACS  (Sorting  of  cells  will  take  place  within 
hours  after  collection)  (2-12) 

Tissue  acquired  via  bronchial  brushing  Is  dissociated,  dead  cells  and  red  blood  cells  are  excluded 
via  FACS,  and  live  cells  are  sorted  Into  96-well  plates  and  frozen.  This  process  needs  to  be 
completed  Immediately  after  sample  collection.  Therefore,  we  have  sorted  all  samples  collected 
above  and  completion  Is  dependent  on  sample  collection. 

Completion  Percent:  33% 

■  Subtask  4:  mRNA  isolation  and  library  preparation  (2-1 3) 

An  established  technique  (CEL-Seq)  for  preparing  single  cell  RNA-Seq  libraries  was  adapted  for 
this  project  and  modified  to  Increase  sample  multiplexing  capacities  and  correct  for  experimental 
amplification  biases.  To  date,  we  have  produced  high  quality  single  cell  data  using  this  protocol 
and  In  the  coming  year  we  will  process  the  samples  collected.  Completion  Percent:  33% 

■  Subtask  5:  Sequencing  of  samples  on  llumina  HiSeq  2500  (1 2-1 4) 

Single  cell  RNA  libraries  are  massively  multiplexed  and  paired-end  sequencing  Is  performed 
using  the  lllumlna  HISeq  2500.  Currently  we  have  only  sequenced  cells  from  our  pilot 
experiments.  Completion  Percent:  15% 

■  Subtask  6:  De-multiplex  samples,  preprocess,  align,  and  analyze  data  quality  (12-18) 

A  computational  pipeline  developed  In  collaboration  with  the  Yana!  Lab 
(http://yanallab. technion. ac.ll/)  was  used  to  preprocess  and  align  reads  generated  via  the  CEL- 
Seq  methodology.  Additional  metrics  have  been  Incorporated  for  the  purposes  of  quality  control 
and  determination  of  sample  cell  type.  The  pipeline  developed  will  be  used  to  process  the  data 
as  It  Is  generated.  Percent  complete:  25% 


■  Milestone(s)  Achieved:  Generation  of  high  quality  single  cell  RNA  sequencing  data  from  30 
subjects 

As  stated  above,  due  to  delayed  IRB  approval  we  didn’t  start  collect  samples  when  we  expected. 
In  order  to  quickly  collect  the  numbers  of  samples  proposed  we  are  adding  another  collection  site. 
Once  the  samples  are  collected,  all  the  protocols  are  now  In  place  to  process  them.  We  expect 
completion  of  this  milestone  within  the  next  year. 


Major  Task  2:  Determine  the  cell  type(s)  responsible  for  the  aberrant  gene  expression  in  the  airway  lung 
cancer  diagnostic  biomarker.  Completion  Percent:  0% 

■  Subtask  1:  Summarize  sequencing  data  into  counts  per  gene  (12-18) 

■  Subtask  2:  Classify  gene  expression  as  signal  or  noise  based  on  a  mixture  model  (12-1 8) 

■  Subtask  3:  Determine  cell  types  of  origin  for  lung  cancer  biomarker  genes  (12-20) 

■  Subtask  4:  Identify  cell  type  dependent  lung  cancer  associated  differential  expression  using  a 
linear  modeling  and  ANOVA  strategy  (18-22) 
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■  Subtask  5;  Choose  candidates  for  validation  (22) 

■  Milestone(s)  Achieved:  Identification  of  which  cell  types  in  the  airway  epithelium  express  lung- 
cancer  specific  gene  expression  alterations 

•  Specific  Aim  2:  Identify  unique  cell  populations  in  the  airway  of  smokers  that  are  associated  with  lung  cancer. 

Completion  Percent:  0% 

0  Major  Task  1 :  Molecularly  identify  subpopulations  of  cells  irrespective  of  cell  type  and  determine  if  these 
cells  are  more  or  less  abundant  in  the  airways  of  patients  with  lung  cancer 

■  Subtask  1 :  Identification  of  novel  subpopulations  of  cells  using  both  class  discovery  and  pathway 
prediction  approaches  (18-22) 

■  Subtask  2:  Identify  subpopulations  associated  with  lung  cancer 

■  Subtask  3:  Choose  candidates  for  validation  (18-22) 

■  Milestone(s)  Achieved:  Identification  of  novel  subpopulations  of  airway  epithelial  cells  that  are 
associated  with  lung  cancer  (22) 

0  Major  Task  2:  Validate  lung  cancer  associated  genes  from  specific  cell  types  or  from  novel 
subpopulations  in  bronchial  epithelial  cells  from  independent  subjects  (n=1 0)  using  FISH 

■  Subtask  1 :  Collect  airway  brushes  from  1 0  subject  for  validation  (1 3-22) 

■  Subtask  2:  RNA-FISH  will  be  used  to  validate  5  candidate  genes  in  conjunction  with  4  known 
epithelial  marker  genes  (22-24) 

■  Milestone(s)  Achieved:  Validation  of  novel  lung  cancer-associated  gene  candidates  in  specific 
populations  of  epithelial  cells 


What  was  accomplished  under  these  goals? 

The  major  objective  is  to  conduct  single  cell  RNA  sequencing  on  cells  collected 
from  bronchial  brushings  from  current  and  former  smokers  with  and  without 
lung  cancer.  In  order  to  attain  this  goal  several  important  experimental 
protocols  needed  to  be  developed  and  validated.  We  have  developed  an 
unbiased  methodology  for  sorting  the  cells  from  bronchial  brushings  into  96- 
well  plates,  a  library  preparation  protocol  that  will  produce  high  quality  data, 
and  an  analysis  pipeline  for  processing  the  data  (Figure  1).  In  order  to 
accomplish  these  tasks  we  have  and  are  continuing  to  conduct  sequencing 
experiments  to  validate  that  our  methods  are  working  prior  to  processing 
brushings  from  smokers  with  and  without  lung  cancer. 

Sample  Collection 

Subjects  recruited  are  current  and  former  smokers  undergoing  flexible 
bronchoscopy  for  clinical  suspicion  of  lung  cancer  at  BUMC.  For  each 
consented  subject,  we  collect  data  regarding  their  age,  gender,  race,  and  a 
detailed  smoking  history  (Table  1). 

Additional  samples  were  collected  from  healthy  current  and  former  smoker 
volunteers  recruited  for  a  project  entitled  “Airway  Epithelium  Profiling  for 
Evaluation  of  E-cigarettes  &  Tobacco  Products”  (Table  1).  The  subjects  in  this 
study  were  recruited  via  advertisements,  a  two-part  screening  questionnaire 
was  administered  via  telephone  to  determine  eligibility,  and  eligible  subjects 
completed  two  study  visits,  the  second  of  which 
included  a  bronchoscopy.  These  healthy  donors 
will  be  treated  as  controls  and  data  generated  for 
the  “Airway  Epithelium  Profiling  for  Evaluation  of  E- 
cigarettes  &  Tobacco  Products”  project  will  be  used 
in  conjunction  with  data  generated  from  samples 
collected  for  this  project. 


Collection  of  bronchial  brushings  is  done  using  the 
same  technique  as  in  our  prior  studies®’®.  Following 
topical  anesthesia  of  the  upper  airway  using  2% 
lidocaine,  a  bronchoscope  is  introduced  to  the  right 
mainstem  bronchus  and  epithelial  cells  are 
obtained  using  an  endoscopic  cytobrush. 


Collect  tissue  via 
bronschoscopic  brushing 

Tissue  dissociation 
Staining  for  FACS 

Sort  single  live  airway  cells 

I 

CEL-Seq  library  prep 

I 

HiSeq  paired  end  sequencing 

I 

Data  pre-processing  and  QC 

I 

Gene  expression  analysis 


Figure  1.  Flow  diagram  depicting  steps 
involved  in  generation  and  analysis  of  single 
cell  RNA  sequencing  data  from  bronchial 
airway  cells. 


Clinical  Suspicion 
of  Lung  Cancer 

Flealthy  Volunteer 

Cancer  Diagnosis 

5  Cancer,  3  No  Cancer,  1  TBD 

21  No  Cancer 

Smoking  Status 

5  Current,  4  Former 

11  Current,  10  Former 

Age 

65  (5) 

44  (11) 

Pack  years  (SD) 

38  (34) 

11  (16) 

Sex 

3  Male,  6  Female 

10  Male,  11  Female 

Table  1.  Demographic  and  clinical  information  from  bronchial  brushing  donors. 
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Cell  sorting  by  FACS 

Tissue  obtained  from  bronchiai  brushings  is  treated  with  0.25%  Trypsin/EDTA  for  epitheiiai  sheet  dissociation  and  ceiis 
are  sorted  using  a  BD  FACSAria  II.  Gating  based  on  forward  scatter  height  vs.  forward  scatter  area  (FSC-H  vs.  FSC-A)  is 
appiied  to  sort  oniy  singiet  events  (singie  ceiis).  Staining  for  GYPA  (CD235a)  is  used  to  exciude  aii  red  bioods  ceiis. 
Hoechst  33342  is  used  to  stain  the  DNA  of  aii  ceiis,  whereas  PI  is  used  to  specificaiiy  stain  the  DNA  of  dead  ceiis  with 
compromised  membranes.  For  each  donor,  singie  Hoechst  33342'"  PI'  CD235a’  ceiis  are  sorted  into  five  96-weii  PCR 
piates  (480  ceiis),  frozen  on  dry  ice  and  stored  at  -80  °C  untii  preparation  for  sequencing  (Figure  2). 
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Figure  2.  Sorting  of  bronchiai  airway  ceiis.  FACS  is  used  to  (A)  isolate  singlet  events  based  on  forward  scatter 

height  vs.  forward  scatter  area  {FSH-H  vs.  FSH-A).  (B)  Dead  cells  (PIT  and  red  blood  cells  (GYPA*)  are  stained  and 

excluded,  and  (C)  live  cells  (Floechst  33342*)  are  sorted. 

Development  of  a  single  cell  mRNA  sequencing  protocol 

Massiveiy  paraiiei  singie  ceii  RNA-sequencing  of  human  bronchiai  airway  ceiis  is  being  performed  using  the  CFL-Seq 
RNA  iibrary  preparation  protocoi  that  has  been  modified  for  increased  high  throughput  capacities”.  Frozen  96-weii  PCR 
piates  containing  sorted  ceiis  are  thawed  on  ice  and  nucieic  acids  are  purified  from  ceii  iysates  using  RNACiean  XP 
beads.  RNA  is  reverse  transcribed  (Ambion  AM1751)  using  primers  composed  of  an  anchored  poiy(dT),  the  5’  liiumina 
adaptor  sequence,  a  weii-specific  barcode,  a  random  sequence,  which  serves  as  a  transcript-specific  unique  moiecuiar 
identifier  (UMI)^^  and  a  T7  RNA  poiymerase  promoter.  Sampies  are  additionaiiy  suppiemented  with  FRCC  RNA  Spike-In 
mix  (Ambion)  for  quaiity  controi^^.  cDNA  generated  from  each  of  the  96  ceiis  per  piate  is  uniqueiy  barcoded  and  therefore 
can  be  pooied  for  second  strand  synthesis  (Ambion  AM1751)  and  ampiification  by  in  vitro  transcription  (Ambion  AM1751). 
Ampiified  RNA  is  then  chemicaiiy  fragmented  (NFB  E61 50)  and  iigated  to  the  liiumina  RNA  3'  adapter  (liiumina  RS-200- 
0012).  Sampies  are  again  reverse  transcribed  and  ampiified  using  indexed  liiumina  RNA  PCR  primers  (liiumina  RS-200- 
0012).  This  barcoding  strategy  aiiows  for  the  ioading  of  iibraries  derived  up  to  4608  ceiis  (96  weii-specific  barcodes  x  48 
piate-specific  indices)  onto  a  singie  fiow  ceii  iane  for  paired-end  sequencing  using  the  liiumina  HiSeq  2500. 

Development  of  a  computational  pipeline  to  process  sequencing  data 

In  collaboration  with  the  Yanai  lab  (http://vanailab.technion.ac.il/),  we  utilized  a  computational  pipeline  to  preprocess  and 
align  reads  generated  from  the  CFL-Seq  protocol.  Briefly,  reads  were  demultiplexed  into  “plate”  level  FASTQ  files  using 
liiumina  software.  Each  plate  FASTQ  file  was  further  demultiplexed  into  cell  specific  FASTQ  files  using  custom  scripts. 

The  UMI  barcode  was  also  removed  from  each  read  and  stored  in  the  read  header.  Reads  were  aligned  to  the  human 
genome  using  Bowtie2.  Gene  level  counts  were  derived  using  a  modified  version  of  the  HTSeq  python  library,  which  only 
counts  multiple  reads  that  are  aligned  to  the  same  position  and  have  the  same  UMI  once,  thus  eliminating  PCR 
amplification  bias.  Quality  of  each  cell  was  assessed  by  examining  the  total  number  of  reads  aligned,  the  total  number  of 
genes  detected,  the  number  of  FRCC  spike  in  transcripts  detected,  and  known  airway  cell  type  markers. 

Pilot  Sequencing  Experiment 

A  pilot  experiment  was  designed  to  determine  if  the  described  single  cell  RNA  sequencing  approach  can  be  used  to 
effectively  detect  a  response  to  cigarette  smoke  exposure  at  single  cell  resolution  and  to  assess  whether  or  not  this 
method  generates  significant  cell-to-cell  technical  variation.  Bronchial  epithelial  cells  from  1  healthy  former  smoker  and  1 
healthy  current  smoker  (n=24  cells  per  donor)  were  profiled  (48  total  cells)  via  paired  end  sequencing  using  the  liiumina 
MiSeq.  Epithelial  cells  were  sorted  using  the  FACSAria  II,  according  to  expression  of  the  established  airway  epithelial 
surface  marker  00166^"'.  Expression  of  major  epithelial  subpopulation  markers  was  detected  in  subsets  of  cells  profiled 
(KRT5,  MUC5AC,  SCGB1 A1 ,  FQXJ1 ).  A  surface  marker  was  used  for  cell  sorting  in  this  pilot  experiment  to  specifically 
generate  RNA  sequencing  data  from  airway  epithelial  cells,  but  all  future  experiments  involve  unbiased  sorting  of  live  cells 
to  profile  all  cellular  subpopulations  present  within  bronchial  brushings. 


7 


To  assess  technical 
variation  of  the 
sequencing  protocol,  the 
ERCC  RNA  Spike-In  Mix 
(Ambion  4456740)  is 
added  to  RNA  from  each 
individual  cell  prior  to 
beginning  CEL-Seq 
library  preparation.  The 
ERCC  RNA  Spike-In  Mix 
is  a  set  of  RNA  controls 
added  at  a  range  of 
known  concentrations. 

The  mean  detection  of 
ERCC  transcripts  across 
all  48  cells  strongly 
correlates  (Pearson’s  r  = 

0.92)  with  known 
concentrations  (Fig.  3A). 

Results  indicate  that  this 
sequencing  approach 
can  be  used  to  generate 
data  with  negligible  cell- 
to-cell  technical  variation. 

Negative  binomial  generalized  linear  modeling  was  used  to  identify  genes  whose  expression  changed  with  respect  to 
smoking  status  (FDR  q  <  0.05).  Differentially  expressed  genes  were  intersected  with  a  set  of  genes  previously  reported  to 
change  in  the  airways  of  current  and  former  smokers  (Fig.  3B)^  All  overlapping  genes  (n=13)  directionally  changed  in 
concordance  with  previously  reported  findings.  Additional  donors  and  increased  numbers  of  cells  per  donor  are  necessary 
to  properly  assess  the  single  cell  resolution  response  to  smoking,  but  these  preliminary  findings  illustrate  that  it  is  likely 
that  single  cell  RNA  sequencing  will  allow  us  to  understand  the  cell  type-specific  contributions  to  the  lung  cancer  field  of 
injury. 

What  opportunities  for  training  and  professionai  development  has  the  project  provided? 

Grant  Duclos,  a  graduate  student,  responsible  for  sorting  the  cells  and  preparing  the  libraries  as  part  of  the  project  has 
had  the  opportunity  to  study  under  Dr.  Itai  Yanai,  Associate  Professor  at  Technion  -  Israel  Institute  of  Technology.  Dr. 
Yanai  is  on  sabbatical  at  the  Broad  Institute  and  has  extensive  experience  in  single  cell  sequencing  and  under  his 
mentorship  Grant  has  been  able  to  effectively  troubleshoot  our  protocols. 

How  were  the  results  disseminated  to  communities  of  interest? 

Nothing  to  Report. 

What  do  you  plan  to  do  during  the  next  reporting  period  to  accomplish  the  goals? 

In  order  to  accomplish  the  goals  in  the  Statement  of  Work,  over  the  next  year  we  plan  to  complete  our  sample  collection. 
We  will  accomplish  this  by  continuing  to  collect  samples  from  subjects  undergoing  bronchoscopy  for  suspicion  of  lung 
cancer  at  Boston  University  Medical  Center.  We  also  plan  to  collect  additional  samples  through  collaboration  with  Dr. 
Robert  Browning,  an  interventional  pulmonologist  at  Walter  Reed  National  Military  Medical  Center.  Grant  Duclos  from  our 
group  will  travel  to  Maryland  to  train  Dr.  Browning’s  group  to  collect,  sort,  and  freeze  cells  obtained  at  bronchoscopies  that 
will  be  shipped  to  BU  for  processing.  Currently,  we  are  in  the  process  of  amending  our  IRB  to  include  this  additional  site. 
Additionally,  we  have  shown  above  that  we  have  successfully  created  protocols/methods  for  collecting,  sorting, 
processing,  and  analyzing  the  data;  therefore,  within  the  next  year  we  plan  to  generate  data  on  approximately  200  cells 
from  about  30  smokers  with  and  without  lung  cancer.  We  will  analyze  the  data  to  discover  cell  type  dependent  lung 
cancer-associated  gene  expression  alterations  in  known  and  unique  cell  populations.  In  addition,  we  will  continue  to 
collect  samples  for  RNA-FISH  validation. 

IMPACT: 

What  was  the  impact  on  the  development  of  the  principal  discipline(s)  of  the  project? 

To  date,  human  bronchial  epithelial  cells  obtained  via  bronchoscopy  have  not  been  sorted  and  processed  for  single  cell 
RNA  sequencing.  The  cells  need  to  be  immediately  taken  from  the  bronchoscopy  suite  and  processed  for  FACS  sorting. 


rV,  rrVTrT ^  i\  h  ffTpn 


Log^Q  ERCC  input  concentration 


J  Former  smoker  cells 
I  Current  smoker  cells 


Expression  level  z-score 


Figure  3.  RNA-Seq  profiling  of  single  bronchial  epithelial  cells  from  healthy  former  and  current  smokers. 

Single  bronchial  epithelial  cells  from  1  former  smoker  and  1  current  smoker  (n=24  cells  each)  were  profiled  using  the 
CEL-Seq  protocol.  (A)  Mean  detection  levels  of  ERCC  RNA  Spike-In  transcripts  across  all  48  cells  strongly  correlated 
(Pearson’s  r  =  0.92)  with  known  concentrations  of  each  transcript  added  to  each  sample.  (B)  Differentially  expressed 
genes  (FDR  q  <  0.05)  associated  with  smoking  status  overlapped  with  and  are  changing  concordantly  with  13  genes 
previously  reported  to  change  in  the  airways  of  current  smokers  [ref  =  Beane  et  al.  2007]. 
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FACS  sorting  of  the  cells  is  accomplished  as  stated  above  and  the  sorted  cells  are  frozen  in  96-well  plates.  The  CEL-Seq 
RNA  preparation  protocol  has  been  modified  to  provide  the  ability  to  process  hundreds  if  cells  from  several  subjects.  All 
of  the  methodology  and  protocols  developed  as  part  of  this  project  will  directly  benefit  other  groups  that  are  attempting  to 
examine  single  cell  transcriptomics  in  human  clinical  samples.  In  addition,  all  microarray  or  RNA  sequencing  studies  to 
date  on  human  airway  epithelial  cells  have  been  conducted  using  a  bulk  population  of  cells.  Expression  values  represent 
the  mean  behavior  of  all  the  cells  profiled  in  a  given  sample  and  do  not  capture  cell-to-cell  gene  expression  variation. 

Using  a  single  cell  approach,  we  will  be  able  to  characterize  the  cell  types  (both  known  and  novel)  that  are  present  within 
the  subjects  and  identify  the  cell  types  responsible  for  the  lung  cancer-specific  airway  field  of  injury.  Using  our  sorting 
technique,  we  will  also  be  profiling  immune  cells  present  when  the  brushing  was  obtained.  Profiling  of  both  epithelial  and 
immune  cells  will  allow  us  to  look  at  the  interaction  between  the  cell  types  and  begin  to  characterize  the  impact  of  the 
immune  system  in  maintaining  a  healthy  epithelium.  The  study  will  contribute  to  our  understanding  of  the  human  airway 
epithelium  and  the  changes  that  occur  in  response  to  smoke  exposure  and  the  acquisition  of  lung  cancer. 

What  was  the  impact  on  other  disciplines? 

The  project  is  multidisciplinary  and  is  likely  to  impact  the  study  of  lung  cancer  and  epithelial  cell  biology  as  well  as 
contribute  to  molecular  biology  and  bioinformatics  methods.  An  example  of  a  broad  impact  would  be  if  the  study  suggests 
new  markers  to  study  rare  epithelial  cell  types  that  appear  to  be  important  in  the  process  of  lung  carcinogenesis.  Another 
example  would  be  if  bioinformatics  techniques  developed  during  the  analysis  of  the  data  for  the  project  could  be  applied  to 
similar  datasets. 

What  was  the  impact  on  technology  transfer? 

Nothing  to  report. 

What  was  the  impact  on  society  beyond  science  and  technology? 

We  have  developed  a  gene  expression  based  biomarker  for  improving  lung  cancer  diagnosis  known  as  PERCEPTA™ 
that  has  been  commercialized  by  Vercyte  (http://www.veracvte.com).  The  test  helps  identify  patients  at  low  risk  for  having 
lung  cancer  after  a  non-diagnostic  bronchoscopy  ordered  as  a  result  of  CT  scan  abnormalities.  The  findings  in  this  project 
may  help  establish  a  more  sensitive  diagnostic  test  that  will  impact  the  clinical  management  of  high  lung  risk  patients. 

CHANGES/PROBLEMS: 

Changes  in  approach  and  reasons  for  change 

We  have  previously  gained  approval  to  switch  single  cell  sequencing  protocols  from  SMART-Seq  to  CEL-Seq  in  order  to 
increase  the  total  number  of  cells  that  will  be  profiled  per  patient.  Using  this  protocol  allowed  us  to  switch  to  a  new  sorting 
procedure  that  captures  all  epithelial,  immune,  and  potentially  unknown  cell  types  in  an  unbiased  fashion.  We  have  also 
expanded  sample  collection  to  include  bronchial  brushings  collected  from  healthy  current  and  former  smoker  volunteers 
recruited  for  a  project  entitled  “Airway  Epithelium  Profiling  for  Evaluation  of  E-cigarettes  &  Tobacco  Products”.  These 
healthy  donors  will  be  treated  as  controls  and  data  generated  for  the  “Airway  Epithelium  Profiling  for  Evaluation  of  E- 
cigarettes  &  Tobacco  Products”  project  will  be  used  in  conjunction  with  data  generated  from  samples  collected  for  this 
project.  Determining  which  cell  types  the  lung  cancer  biomarker  genes  are  expressed  in  these  additional  control  groups 
will  increase  our  understanding  of  the  nature  of  the  “field  of  injury”  and  why  these  genes  are  differentially  expressed  in 
lung  cancer  patients  only. 

Actual  or  anticipated  problems  or  delays  and  actions  or  plans  to  resolve  them 

With  CEL-Seq,  we  are  detecting  -2000  genes  per  cell  on  average.  It  is  possible  that  this  coverage  will  not  be  enough  to 
distinguish  all  subpopulations  of  cells  within  the  airway  brushings.  If  this  occurs,  we  will  either  perform  deeper  sequencing 
of  the  CEL-Seq  libraries  to  detect  more  rare  transcripts  or  perform  SMART-Seq  on  an  additional  set  of  cells  which  is 
reported  to  detect  -5000  genes  on  average^®. 

In  addition,  we  have  had  delays  in  sample  collection  because  of  the  IRB  approval  process.  We  will  continue  collecting 
samples  from  Boston  University  Medical  Center  and  will  also  collect  brushes  in  collaboration  with  Dr.  Robert  Browning’s 
group  at  Walter  Reed  Medical  Center  as  explained  above. 

Changes  that  had  a  significant  impact  on  expenditures 

Nothing  to  Report. 

Significant  changes  in  use  or  care  of  human  subjects,  vertebrate  animais,  biohazards,  and/or  seiect  agents 

Nothing  to  report. 

Significant  changes  in  use  or  care  of  human  subjects 

Nothing  to  report. 


Significant  changes  in  use  or  care  of  vertebrate  animals. 

Nothing  to  report. 


Significant  changes  in  use  of  biohazards  and/or  select  agents 

Nothing  to  report. 

PRODUCTS: 

Publications,  conference  papers,  and  presentations 
Nothing  to  report. 

Journal  publications. 

Nothing  to  report. 

Books  or  other  non-periodical,  one-time  publications. 

Nothing  to  report. 

Other  publications,  conference  papers,  and  presentations. 

Nothing  to  report. 

Website(s)  or  other  Internet  site(s) 

Nothing  to  report. 

Technologies  or  techniques 

The  techniques,  inciuding  ceii  sorting,  iibrary  preparation,  and  anaiysis  methods  deveioped  in  this  project  wiii  be  shared 
through  pubiication  of  a  manuscript  reporting  the  findings  of  singie  ceii  sequencing  experiments  profiiing  smokers  with  the 
without  iung  cancer. 

Inventions,  patent  applications,  and/or  licenses 

Nothing  to  report. 

Other  Products 

Nothing  to  report. 

PARTICIPANTS  &  OTHER  COLLABORATING  ORGANIZATIONS 
What  individuals  have  worked  on  the  project? 


Name; 

Jennifer  Beane-Ebel 

Project  Roie; 

Principal  Investigator 

Researcher  Identifier  (e.g. 
ORCID  ID): 

Nearest  person  month 
worked; 

3 

Contribution  to  Project: 

Overseeing  all  aspects  of  the  project  including  sample  collection, 
experimental  design,  and  data  analysis 

Funding  Support: 

DOD  IDA,  LUNGevity  Foundation,  NIH/NIAID,  Industry  award. 
Internal  funds 

Name; 

Joshua  Campbell 

Project  Role; 

Postdoctoral  Fellow 

Researcher  Identifier  (e.g. 

10 


ORCID  ID): 


Nearest  person  month 
worked: 

5 

Contribution  to  Project: 

Overseeing  aii  aspects  of  the  project  with  the  PI  and  leading  the 
computational  analysis  of  the  data 

Funding  Support: 

DOD  IDA,  NIH/NHLBI,  Industry  funds 

Name: 

Martine  Dumas 

Project  Role: 

Study  Coordinator 

Researcher  Identifier  (e.g. 
ORCID  ID): 

Nearest  person  month 
worked: 

1 

Contribution  to  Project: 

Consenting  patients  and  sample  collection 

Funding  Support: 

DOD  IDA,  NIH/NCI,  Industry  funds 

Name: 

Grant  Duclos 

Project  Role: 

Graduate  Student 

Researcher  Identifier  (e.g. 
ORCID  ID): 

Nearest  person  month 
worked: 

6 

Contribution  to  Project: 

Collecting  samples,  sorting  cells,  and  preparing  RNA  sequencing 
libraries  from  the  cells,  and  data  analysis 

Funding  Support: 

DOD  IDA,  NIH/NHLBI 

Name: 

Yaron  Gesthalter 

Project  Role: 

Researcher  Identifier  (e.g. 
ORCID  ID): 

Nearest  person  month 
worked: 

NIH/NI 

Contribution  to  Project: 

Consenting  patients  and  collection  of  samples 

Funding  Support: 

NIH/NHLBI,  Internal  funds 

Has  there  been  a  change  in  the  active  other  support  of  the  PD/Pi(s)  or  senior/key  personnei  since  the  iast 
reporting  period? 

Nothing  to  Report 
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What  other  organizations  were  involved  as  partners? 

Organization  Name;  Broad  Institute 

Location  of  Organization:  Boston,  MA 

Partner’s  contribution  to  the  project:  Coiiaboration.  We  have  been  coiiaborating  with  Itai  Yanai,  Associate  Professor  at 

Technion  -  Israei  Institute  of  Technoiogy  during  his  sabbaticai  at  the  Broad  Institute  to  heip  deveiop  our  singie  ceii 

sequencing  iibrary  preparation  protocoi. 

SPECIAL  REPORTING  REQUIREMENTS 

None. 
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